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Abstract 



We report a glueball Regge trajectory emerging from diagonalizing a confining 
Coulomb gauge Hamiltonian for constituent gluons. Using a BCS vacuum 
ansatz and gap equation, the dressed gluons acquire a mass, of order 800 MeV, 
providing the quasiparticle degrees of freedom for a TDA glueball formulation. 
The TDA eigenstates for two constituent gluons have orbital, L, excitations 
with a characteristic energy of 400 MeV revealing a clear Regge trajectory 
for J = L + 5, where S is the total (sum) gluon spin. Significantly, the 8 = 2 
glueball spectrum coincides with the Pomeron given by ap(t) = 1.08 + 0.25 t. 
Finally, we also ascertain that lattice data supports our result, yielding an 
average intercept of 1.1 in good agreement with the Pomeron. 



Pacs #: 11.55.Jy; 12.39.Mk; 12.39.Pn 12.40.Yx Keywords: Glueball Regge trajectories, 
Pomeron, QCD, Coulomb gauge Hamiltonian, TDA. 
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A vintage problem in hadronic physics is to fundamentally explain the observed 
monotonously rising, to asymptotically fiat cross section behavior with increasing Man- 
delstam variable s. Pre-QCD cross section theorems involving t channel exchanges mo- 
tivated Pomeranchuk's assertion that vacuum quantum number {J^'-' = O^'^) t channel 
exchanges would describe elastic cross section high energy dependence. Related, it was his- 
torically known that hadronic resonances with specific quantum numbers could be connected 
to a series or "tower" of excited states. From this backdrop the successful Regge picture 
0,^, with poles and trajectories linear in t (mass squared), a{t) = a{0) +bt, emerged which 
effectively described cross sections and resonances and therefore unified both scattering and 
bound state data. Of particular interest in this letter are analyses of high energy diffraction 
data 1^,1^ yielding a J^*" = Regge trajectory with the largest intercept a(0). This 
trajectory is called the Pomeron and from detailed fits is given by 

ap(t) = 1.08 + 0.25 t . (1) 

Note that by definition, the Pomeron entirely governs and correctly reproduces the high 
energy scattering behavior since the cross section scales as s°'^^^~^. The simplicity of this 
picture is very appealing and the Pomeron remains a contemporary tool for understanding 
high energy processes, both experimentally and theoretically [0] . In the discussion below we 
specialize to the soft non-perturbative Pomeron (for an account of the BFKL hard Pomeron 
in perturbative QCD see 

According to Regge theory, the trajectory describing scattering should also be consistent 
with physical hadron states, implying a J^*" = J~^~^ trajectory of mesons with slope near 
0.25 GeV~^. However, all known hadronic resonances (both mesons and baryons) fall on 
Chew-Frautschi []T0[ trajectories with slopes close to 0.9 GeV~^ but different intercepts. This 



is illustrated in Figs. |T| and |^ for several mesons and baryons, respectively, including the p 
whose negative parity, positive G parity tower is well described by ap{t) = .55 + 0.9 GeV~^ t. 
This incompatibility led to the conjecture that the Pomeron corresponded instead to a Regge 
trajectory for gluonic states, and most likely one describing the J"*""*" glueballs. 



Not surprisingly, the gluonic nature of the Pomeron has been widely investigated (see |TT 
for a recent review) with several recent studies linking the Pomeron and glueball trajectories. 
From meson trajectories and Wilson loops, Kaidalov and Simonov [|12[ estimate that the 



glueball intercept is a{fi) ~ .7, while Soloviev |jT3[ computes q;(0) = 1.07 ± 0.03 using 



the quantized elliptic Nambu-Gotto string model. Motivated by the Maldacena conjecture. 



and solving the dilaton wave equation, Brower, Mathur and Tan |14| obtained a glueball 



intercept of a(0) ~ 1.2. Finally, using mass relations between the glueball and meson 



sectors, Brisudova, Burakovsky and Goldman [|1^] obtain a glueball Regge trajectory with 
slope ~ 0.3 ± 0.1 GeV^. 

In this communication we report results from a new constituent glueball approach which 
further strengthens the gluonic Pomeron conjecture. Similar to previous relativistic quark 



approaches p!6|JT7[] yielding a consistent meson Chew-Frautschi Regge slope, our relativistic 



many-body calculations also produce a glueball Regge trajectory that arises naturally from 
the confining linear potential and the scale of angular excitations. In reproducing the lattice 
gauge positive parity glueball spectrum, we obtain a Regge trajectory with intercept near 
and above unity. Most significantly, the trajectory slope is close to the Pomeron value 
{h ^ 0.25 GeV"^) and predominantly independent of model details. Our group has applied 
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many body techniques to develop a unified approach for the hadron spectrum utihzing 



a QCD inspired Hamiltonian. Both the gluon ||18| and meson [|19| , p0| sectors have been 
reahstically described and a hybrid meson calculation has also recently been completed ^1 . 

We briefly summarize our many-body approach (see refs. [p!^-pO| for complete calcula- 
tional details). The starting point is a Coulomb gauge, field theoretical Hamiltonian 

H = TrJ cix(n« • n'^ + ■ B^) - i I rfxrfyp»(x)\/(|x - y|)p«(y) , (2) 

with color charge density p° = f'^^'^A^ ■ and fields 11", B^ = V x A"^. We utilize an 
instantaneous Cornell potential for confinement, V = — — + ar. The potential parameters 
are as = .2, and, from independent lattice results, o" = | 0.18 GeV^. The gluons are then 
dressed by means of a Bardeen, Cooper and Schrieffer (BCS) ansatz for the vacuum ground 
state which, through a variational Hamitonian minimization, generates the gap equation 

+ (3) 

for the gluon self-energy, Ug, containing the dynamical gluon mass. Our potential in mo- 
mentum space is V^(|q — k|) = —-^z^ ~ (q-kY ' ^^^^^ ea.ch. term generates a divergence in 
the gap equation (Coulomb is quadratic, confinement is logarithmic). Since we are focusing 
on confinement we have omitted the Couloub potential component and then used a cut-off 
parameter A = 4 — 5 GeV to regularize the remaining logarithmic divergence. Finally, a 
Tamm-Dancoff (TDA) diagonalization, truncated to the Ip-lh level, is performed for the 
same field theoretical Hamiltonian (no divergences. Coulomb potential is now included) to 
generate the glueball spectrum. We used the same linear potential in our quark sector ap- 
plications which produced reasonable results [|r^,^], including the appropriate meson Regge 
slope as indicated by the solid squares in Fig. |^ for the 1 and 3 p trajectory. It is 
important to stress that our quark sector application properly implemented chiral symmetry. 

At the more technical level, the glueball TDA state, for fixed orbital, L, and total gluon 
spin, 5, is represented by 

l^ir)= E /7^*i£m.(q)<(qX(-q)l^), (4) 
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where the are quasiparticle creation operators acting on the BCS vacuum state \VL). The 
sum is over the color index a and the two transverse spin projections, m = 1, 2. The glueball 
wavefunction has angular momentum composition 

H'^sZ.mM) = {lmam,\Sms){Lm,Sms\Jmj)<pfs''i<l)yri^) ■ (5) 

As mentioned above, the gluon quasiparticle energies are obtained from the nonlinear mass 
gap equation which yields a gluon mass around 800 MeV . Then, with the gap energy, 
the linear TDA equations are diagonalized for the J^^ states of interest, either as a radial 
equation or variationally in a multi-dimensional Monte Carlo calculation for the Hamilto- 
nian matrix elements using the computer code VEGAS. Our work omits renormalization, 
but a recently improved, renormalized glueball calculation produced similar J^*^ = O^"*" 
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results. Using a new, independent code we have recomputed the glueball spectrum and 



confirmed our original results We have also extended the analysis to the maximally ac- 
ceptable Ip-lh excitation range spanning higher angular momentum states, up to S^"*". The 
Ip-lh energy upper bound is about 3 GeV since it is well known that 2p-2h (4 gluon) states 
become important for energies twice the lightest Ip-lh excitation 1.6 GeV ground state 
glueball). A realistic calculation above 3 GeV would entail mixing with 2p-2h states requir- 
ing a complicated (4 gluon) Fock space diagonalization. Further, higher angular momentum 
significantly enhances the interaction spin dependence which is not accurately known. These 
two reasons are also why we only predict states up to J = 3 for the meson spectrum shown 
in Fig. |I[ The positive parity glueball states reliably predicted by our model are compared 
to the previous and most recent lattice measurements |^ in Fig. ^ 

Wavefunction symmetry for two identical bosons requires that L + S must be even. Since 
the Pomeron should carry positive parity and C-parity, L odd is not considered nor three 
gluon states which have odd C-parity. Also, according to Yang's theorem, a combination 
of gluon transversality (Coulomb gauge) and Bose statistics, J = 1 states are forbidden 
for a two gluon system. Further, the J = state is precluded from the Pomeron since 
it necessarily belongs to a lower (daughter) trajectory. Therefore only the J = 2"*""^, with 
L = 0, 5* = 2, and J = 3~^^, with L = 2, = 2, states in Fig. ^ are Pomeron candidates 
and the corresponding glueball Chew-Frautschi diagram is shown in Fig. ^ along with the 
Pomeron trajectory. 

Figure H constitutes the key findings of this study. Notice that our model yields a tra- 
jectory with slope .20 GeV~^ and intercept 1.1, in reasonable agreement with the Pomeron. 
We found little sensitivity of our results to the confining potential parameter as increasing 
the strength by 33 % (i.e. a = |0.18 GeV'^ — >• 0.18 GeV'^), only marginally reduced the 
intercept. Further, and of interest to BFKL PQCD studies, our model predicts the soft 
Pomeron is a correlated two gluon state with total intrinic spin S = 2. 

Equally significant are the lattice data which qualitatively supports our main result. The 
lattice calculations entail some uncertainty and therefore generate a less accurate trajectory 
having slope .15 ± .05 GeV~^. However, the lattice intercept is above unity, although again 
with appreciable error (1.1 ± .5). We therefore conclude that both our model and recent 
lattice data are consistent with a glueball interpretation of the Pomeron. 

Ideally a J = 4++ (and higher) state should also be predicted. This may be possible 
for lattice calculations but, as discussed above, would be quite formidable in our approach. 
Conversely, assuming a glueball Pomeron interpretation we can use the fitted trajectory to 
predict that the 4++ and 5+"'' glueballs should have masses Mj = 2VJ- 1.08 = 3.42 and 
3.96 GeV, respectively. 

Finally, we offer a qualitative explanation for our unified model predictions of meson 
and Pomeron slope trajectories with respective values bgg ^ 0.9 GeV^ and bgg ~ 0.25 GeV^. 
In a simple QCD string model with massless partons separated by a string with length r 
and tension (or energy per unit length) k, the total relativistic mass scales as kr while the 
angular momentum is proportional to |A;r^, which generates a Regge trajectory with slope 
b = 2^.- This string feature is incorporated in our model through the linear confining 
potential (err), yielding effective string tensions, kgg = Cgga and kgg = Ggga, for the meson 
and glueball systems, respectively. Here Ggg = and Ggg = Nc are the respective quark 
and gluon color Casimir operators from the density- density Hamiltonian kernel. Hence for 
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Nc — 3 the meson and glueball trajectory slopes are related by bgg = -(f-hqq — ^hqq. While 
this simple estimate explains much of the meson/glueball slope difference, it is significant 
that our model is able to account for entirely all of the difference, the rest due to BCS 
quasiparticle masses and the Coulomb potential in the TDA. It would be very interesting 
to see if alternative QCD models, particularly the flux tube, can reproduce these results. 

In summary, our relativistic many-body approach provides a glueball Regge trajectory 
similar to the Pomeron diffraction fit. The model also unifies and describes the quark sector 
by properly implementing chiral symmetry (and breaking) to reproduce the observed meson 
spectra and attending Regge trajectories with the same, predetermined confining potential. 
Recent lattice measurements are also consistent with this picture and support a gluonic 
interpretation of the Pomeron. Reproducing higher states on the Pomeron trajectory 
would confirm this assertion and more sophisticated model calculations are in progress along 
with a call for alternative model predictions. Finally, our results indicate that the Pomeron 
provides a worthwhile guiding constraint which future theoretical and experimental glueball 
studies should utilize. 
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FIG. 1. Meson (p, K* and a) Regge trajectories constructed from recent tabulated data (dark 
circles and error bars, PDG 2000). Boxes are model TDA predictions for the p trajectory. 
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FIG. 2. Baryon (AT and A) Regge trajectories constructed from recent tabulated data (dark 
circles and error bars, PDG 2000). 
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FIG. 3. The TDA (boxes) and lattice gauge (vertical bars showing error) glueball spectra. The 
lattice results, from left to right, are from refs. p3|-p^, respectively. 
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FIG. 4. Comparison of TDA (boxes) and lattice (horizontal bars) J"*"*" glueballs with the 
Pomeron (solid line). 
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